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Abstract 

The one-dimensional Hubbard model with different on-site interactions 
is investigated by renormalization group technique. In the case of a 1/4- 
filled band the dynamical nonequivalence of sites leads to the appearance of 
Umklapp processes in the system and to the dynamical generation of a gap in 
the charge excitation spectrum for Ua 7^ Ub, Ua > or Ub > 0. The ground- 
state phase diagram is obtained in the limit of second order renormalization. 
Depending on the sign and relative values of the bare coupling constants, 
there is a gap in the spin or charge excitation spectrum and the model system 
tends to superconducting or antiferromagnetic order at T = 0, with doubled 
period. The role of interaction between particles on nearest and next-nearest 
neighbor sites is also considered. 

I. INTRODUCTION 

The remarkable discovery of high-Tc superconductivity |l| in Cu — O compounds has 
raised the interest in the physics of correlated electrons in low dimensions. Although the 
main interest is centered around the physics of two-dimensional highly correlated electron 
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systems, the one-dimensional analogues of models related to high temperature superconduc- 
tivity are very popular due to the conjecture that properties of the ID and 2D variants 
of certain models have common aspects. 

Emery [^] proposed a one dimensional version of two-band copper-oxide model In 
Emery's model two essential features distinguishing the Cu and O sites in real systems are 
incorporated: 1. the chemical aspect - the hole has different on-site energies on copper and 
oxygen sites and 2. the dynamical aspect - there are different on-site interactions between 
two holes on copper sites {Ua) and on oxygen {Ub)- Detailed analysis of the two-band 1/4- 
filled Peierls-Hubbard model in the large U limit is given in and in the case of large 

dimerization in [Q. To investigate the influence of the dynamical nonequivalence of copper 
and oxygen sites Japaridze etal. considered a rather formal ID version of the Cu — O 
chain, namely the one-band Hubbard model with different on-site interaction on even and 
odd sites. This model was investigated in the framework of standard weak coupling approach 
using a bosonization technique. 

The another method by which insight can be gained into the problem of interacting 
fermions in ID is the renormalization group approach P|-pil|. In this paper we consider the 



quarter-filled band case of the one-band Hubbard model with alternating interactions using 
the renormalization group technique. 

As will be shown in 2, the continuum Hamiltonian of the model has a general form, de- 
scribing the one-dimensional system of interacting fermions. The Hamiltonian of the same 
form has been considered by Kimura [§] using the conventional field-theoretical renormaliza- 
tion group. Solyom performed the second order scaling procedure, using the bandwidth 
cutoff, to more general case, when the forward scattering process between particles from 
the same branch was presented. According to Solyom, taking into account this process sub- 
stantially changes results. Unlike 0, Rezayi etal. have shown, that if one modifies 
the term corresponding backward scattering and Umklapp process, replacing them by ef- 
fective spin and charge carrying interactions, one gets the theory with momentum transfer 
cutoff. In this case, the effect of the forward scattering process between particles from the 



same branch is absorbed to a Fermi velocity renormalization and couplings, which couple 
the spin-density degrees of freedom scale independently of the couplings, which couple the 
charge- density degrees of freedom. 

In our calculation, we shall neglect the forward scattering process between particles from 
the same branch and assuming that the forces are short range, use the bandwidth cutoff. 
In this case, the Lie equations for the coupling constants decouple into two sets and are 
analogous to those obtained by using the momentum transfer cutoff [0]. As will be shown 
in 4 the scattering processes responsible to a Fermi velocity renormalization can only lead to 
the difference between the degrees of the divergency of CDW and SDW response functions. 

II. THE MODEL 

The Hamiltonian describing the model has the following form: 

H=-tJ2 <a(Cn+l,a + C„_l,„) + ^Y^P + {~^TV]pn,aPn,-a+ (l) 
n,a n,a 

+ 51 [^lPn,aPn+l,l3 + V2pn,aPn+2,l3] 
n,a,f3 

where c^a{^n,a) is a creation (annihilation) operator for a particle of spin a in the Wannier 
state locahzed at the nth lattice site, pn^a = c^a^n,a and 

t/ = i(f/, + f/,), V=]^{Ua~U,) (2) 

where UaiUh) describes on-site interaction on even (odd) lattice sites and Vi and V2 describe 
the interactions between particles on nearest and next-nearest neighbor sites, respectively. 

The continuum limit of the Hamiltonian (1) can be taken by approximating the spectrum 
—2tcosk by linear spectra around ^kp and using the correspondence: 

Cn,a exp{ikpn)ipi^a{x) + exp{—ikFn)tp2,a{x) (3) 

where the fermion field ipi{2),a{x) describes the particles around Fermi point +{—)kp. 

The information about the lattice structure of the initial model (1) is contained in the os- 
cillating factors exp{±ikFn) and manifests itself in the continuum limit Hamiltonian via the 
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selection of the scattering processes relevant for the physics of the system at large distances. 
Except for the cases of quarter- and half-filled bands the dynamical nonequivalence of sites 
is irrelevant (in the sense used in connection with the renormalization group) and the con- 
tinuum limit Hamiltonian has the structure of the well known Luther-Emery backscattering 
model with parameters depending on U In what follows we will consider the case 
of quarter-filled band. In this case the dynamical nonequivalence of sites becomes important 
and the continuum limit Hamiltonian has the form of backscattering model supplemented 
by Umklapp processes with bare coupling constant proportional to V. 

H = -iVpYl / dxllpt^^d^tpi^a - i't,a9xi'2,a] + (4) 



+ [92\\Sa,f3 + 92±Sa,-p]'ipla'^lp^2,l3^l,a + -jSa,-f3{ijt,a'^t,f3^2,p^2,a + h.C.)| 

The Hamiltonian (3) has a general form describing the interacting fermion system in ID. The 
terms with coupling constants gi\\, gi± and g2\\, g2± correspond to the backward and forward 
scattering processes,respectively and the terms with g^ - to the Umklapp process (see Fig. 
1). In the (3) the processes with coupling constants gi\\ and g2\\ are indistinguishable and 
we can choose g2\\ = g2±, without loss of generality. This leads to the following values of the 
spin independent bear coupling constants: 

g^ = U -2V2 g2 = U + 2Vi + 2V2 g3 = V 

In obtaining (3), the terms corresponding to forward scattering processes, in which both 
incoming particles are from the same branch, were omitted. For the following discussion the 
parameter gp = gi — 2g2 is also useful. The coupling constants gs and {gp,gz) describe the 
spin and charge degrees of freedom respectively ll3l . 



4 



III. RENORMALIZATION GROUP TECHNIQUE 



The multiplicative renormalization group has been used in the field theory for a long 
time to eliminate divergences [^. This method has two main aspects. The first one is that 
starting from a perturbative calculation a partial summation is obtained by solving the group 
equations and second one - a set of equivalent problems can be found which are described 
by Hamiltonian of similar form. If there is a model among the equivalent ones which can 
be solved, the solution of the original problem can also be obtained. The renormalization 
group treatment of the problem of interacting fermions in ID was given by Menyhard and 
Solyom |jl2[. Here we will briefiy present the ideas applied to the ID Fermi gas. 



For an unambiguous definition of the model it is necessary to specify the cutoff proce- 
dure which determines the domains of admissible values of the electron momenta. If the 
interactions are short range we have effectively only one limit on momenta and got a theory 
with bandwidth cutoff. In the general case, when the long-range forces are also presented, 
theory needs two cutoffs, one as bandwidth and one on momentum transfer |jl3[ . Assuming 
that the forces are short range in our case, we shall use the bandwidth cutoff. 

At first let us write the auxiliary Green's function d and four-pointed vertex function Fj 
in the following way: 

G{kF, uj) = d{^, gi)Go{kF, uj) (5) 

ra/375('^) = 5'lf l(a;)5a75/35 " 5'2f 2(t^)5a55/37 (6) 

for process with momentum conversation and 

for Umklapp process, where loq is a bandwidth cutoff. 

If the cutoff Eq is changed to and simultaneously the coupling constants gi are trans- 
formed to certain new values 



then the functions d and Fj are multiphcatively transformed to 

d{.^,g'i) = zd{^,gi) (8) 

u^,9';) = z^'u^,9^) (9) 

Moreover, requiring the invariance of quantity giT^d'^, which is an appropriately renormahzed 
vertex [Q, only four of renormalization constants Zi, z[ and z are independent and we have 
A = Zi/Z^. 

The scaling equations for d, Fj and g'^ can be written in a common form as: 

C{^,g',) = zC{^,g,) (10) 



For any quantity obeying the condition (11^) a Lie differential equation of the form: 

(11) 



d d 



can be derived, where x = uj/Eq. 

The right-hand side of Lie equation can be calculated by perturbation theory. If the 
invariant couplings gi{x) = g[ are small for an arbitrary change of the scaling energy the 
quantity determined from the solution of the Lie equation using a few terms of the pertur- 
bation series will present a good approximation in the whole energy range. If, however, the 
invariant couplings become of the order of unity while scaling energy goes towards lower 
energies, the right-side of the Lie equation breaks down and only qualitative results can be 
obtained. 

A. First-order renormalization 

Calculating the right-hand side of the Lie equations in the second-order of the perturba- 
tion theory the following equations for invariant couplings can be obtained 

d 

^g^dsix) = gl{x) (12) 
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d 

^Q^9pix) = glix) (13) 
d 

The first order scaling is equivalent to summing up the leading terms and these equations are 
exactly analogous to those obtained in [|T^ by parquet approximation. Let us first analyze 
the model (1) in the case Vi = V2 = In this case the solutions of Eq. (|l2|) and ( |13| ) with 
boundary conditions gsiS^) = U/t:vf, gp{0) = —U/nvp and (73(0) = V/nvp have the form: 

= ^ (14) 



where t = \nx, 



-U V 
Qoix) = — Qsix) = — (15) 

for U'^ = V'^, while for U'^ > V'^ we have 

g,(.) = -DcotlW(t„-t) ,3(.) = ^J^f|l^ (16) 

where to = l/-Darccoth(f/ /D) and for < they reads as 

= -D.otDit',-t) = (17) 

where = l/Darccot(t//-D). 

These equations show that, at certain values of bare couplings, the poles appear in 
the expressions of the invariant coupling constants. This singular behavior of invariant 
couplings doesn't indicate the phase transition at finite temperature, which is impossible in 
one dimensional system |16|. The existence of the poles in expressions (p!4D-(p!7D indicates 
that at low temperatures the interactions become strong and this approximation is no longer 



valid. In |15| Dzyaloshinsky and Larkin, comparing their result to exact solution obtained 



by Gauden |T^, emphasized, that the poles in the expressions of the invariant couplings 
can indicate the appearance of the gap in the one particle excitation spectrum. Thus we 
can expect that the spin (charge) gap exists when there is a pole in the expression for gg 
((fl'3) fi'p))- Thus the charge gap exists in all sectors except A, while the spin gap exists in 
the sectors A, B, E (see Fig. 2). 



B. Second-order renormalization 



The second-order renormalization corresponds to considering next to leading logarithmic 
corrections and thereby to taking into account fluctuation effects. In the limit of third order 
of perturbation theory for d and Fj, the following Lie equations for invariant couplings can 
be obtained: 

^§^9s{x) = gl{x)[l + ]^gs{x)] (18) 

x^9p{x) = gl{x)[l + ]^gp{x)] (19) 

d 11 

a;^^73(a;) = g,{x)g:i{x)[l + -gp{x)] + -gl{x) (20) 

The solutions of these equations can be obtained only in an implicit form. They are the 
smooth functions of energy and tend to the saturation values, fixed points, at = 0. The 
values of the fixed points can be obtained from zeros of the right-hand sides of the Lie 
equations (|l8D -(pOD using the following arguments: 

d 

a). x—gs{x) < for gs{x) > -2; 
ox 

b). x—gpix) < for > —2 and ] , ^ = c, where c is a constant. 
ox gp[x) + 2 

The values of the fixed points depend on the relation between the bare coupling constants 
and are presented in Table 1. 



C. Response functions and ground-state phase diagram 

The symmetry of the ground-state can be found by investigation of various response 
functions which are expected to be singular. These quantities are the charge- density (CDW) 
and spin-density (SDW) waves with wave-vector k = 2kF, the singlet-superconductor (SS) 



and triplet-superconductor (TS) response functions with wave- vector k = ||T9 |. 
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In the case of 1/4-filled band the periodicity of the CDW and SDW is equal four lattice 
constant and it is possible to introduce two separate sets of order parameters describing the 
fluctuations of the charge and spin density on even (a) and odd (6) sublatticies. In the real 
space and time representation general formula for these functions is 

= -ie{t) < [U{x,t)U+{0,t)] > (21) 

where U{x,t) = J2a[i't,a{^y't)i^2,a{x,t) + ilj2,a{^jt)i^i,oi{x,t)] foT the charge- density response 
function N{x,t); ip^^^{x,t)ip2,i{x,t) + tp2.^{x,t)ipi^i{x,t) for the spin-density response func- 
tion xi^jt) ^]tT(^' ^)^2,i(T)(^' ^) singlet (triplet)-superconductor response function 
As(t)(x,t) 

The Fourier-transformed response functions don't satisfy the criterion of multiplicative 
renormalization, Solyom introduced the auxiliary renormalizable functions defined by 

where ^i{uj) is one of the above considered response functions. In the limit of second order of 
perturbation theory the following Lie equations for these renormalizable response functions 
can be obtained: 



dNa(b){x) 



2gs{x) - g2{x) + {-)gz{x) + F{x) (22) 



(91n(x) 

^^"^ -g2{x) - {+)g,{x) + F{x) (23) 



dXa{b){x) 



d\n.{x) 



^^'^""^ gs{x)+g,{x) + F{x) (24) 



d\n.{x) 

= g2{x)-gs{x) + F{x) (25) 



d^Ax 



where 



d\n.{x) 



F{x) = ^{9six) + 92i^) - 9s{x)g2{x) + i^|(x)) 
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and Na{x) and Nb{x) {Xai^) Xbi^)) describe the CDW (SDW) located on a and b 
sublatticies , respectively and they are different only in a sign of g^, the changing aQb is 
equivalent to change the sign of Umklapp processes. 

In the limit uQO the asymptotic expressions of response functions are obtained as 

$,Hoc$,Hoc(^)"' 

-C/0 

The critical exponents can be obtained by putting the fixed points of invariant couplings 
in the right-hand side of Eqs.(^)- (l25|) and they are presented in Table 2 

According to singularities in the response functions let us summarize the ground-state 
phase diagram. In the sector A the CDW and SS response functions are divergent. The 
dominant singularity is in the singlet-superconductor response. In the sectors B and E only 
CDW response, located on the sublattice with lower on-site interaction, is divergent. There 
is coexistence of CDW and SDW in the sectors C and D. The CDW is still located on the 
sublattice with lower on-site interaction while the SDW is located on the other one. 

In the case of nonalternating chain, = 0, g4 does not renormalize and system has the 
properties of normal metal for U > and singlet-superconductor for U < 0. 

We note ground-state phase diagram obtained here is exactly analogous to that obtained 
by boson representation theory in [0. 



IV. SCALING TO THE EXACTLY SOLUBLE MODELS 

The renormalization group and scaling arguments presented in 3 establish a relationship 
between the original problem and a set of problems in which the coupling constants have 
the somewhat different values. If an exactly soluble model appears among these equivalent 
systems, the physical behaviors of the original model can be predicted by using the scaling 
arguments. 

From Tab.l one can easily see that depending on relation between the bare couplings 
the model can be scaled to Tomonaga |^ or Lutter-Emery (LE) model and the 

ground-state phase diagram can be obtained by using their results ||13| . 
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(i) . In the sector A (see Fig. 2), < 0, gp{l) > \g3{l)\, the spin part of the 
Hamiltonian scales to the LE hne and there is a gap in the spin part of the excitation 
spectrum. The charge part scales to a Tomonaga model and the Umklapp processes have 
no influence. Only the CDW and SS responses can be divergent. The dominant singularity 
is in the latter. 

(ii) . For 5's(l) > 0, gp{l) < \g3{l)\, sectors C, D, the situation is reversed concerning the 
spin and charge parts of the Hamiltonian. The gap is in the charge part and CDW and SDW 
responses show a singular behavior. The latter is more singular. The difference between 



singularity in CDW and SDW responses is caused by Fermi velocity renormalization |T3 
for which scattering processes omitted in (2) is responsible ||10|| . 

(ill). In the sectors B and E, (7^(1) < 0, gp{l) < \g3{l)\, the Hamiltonian scales to 
the LE line, there is a gap in both charge and spin parts and only charge-density response 
function is divergent. 



A. The role of intersite interactions 

Taking into account the effects of intersite interactions leads to the renormalization of 
bare coupling constants Eq.(4). Let us consider the phase diagram for the arbitrary sign 
of the bare coupling constants, not restricting ourselves by realistic case, when Ua, Ub, 
Vi and , 1^2 > 0. 

In the case of extended model the spin and charge gaps exist in the case when U — 2V2 < 
and -{U + 4Vi + 6V2) > \V\, respectively. 

(i) . Now the situation discussed in 4.(i) takes a place when 

[/ - 2^2 < and -{U + AVi + 6V2) > \V\ 

and the ground-state of the system is the same as in sector A. 

(ii) . When the bare couplings constants satisfy the condition: 

f/ - 2F2 > and -{U + AVi + QV2) < \V\ . 
11 



in the extended model, the same response functions as in sectors B and E are divergent, 
(iii). The system has the same ground-state as in sectors C and D when 

[/ - 2^2 < and -(U + Wi + 6V2) < \V\ . 

(iiii). Unhke to the case Vi = V2 = 0, when intersite interactions are presented, the 
Hamiltonian scales to the Tomonaga model when: 

U -2V2>0 and - {U + AV^ + 6V2) > \V\ 

There is no gap in either the charge or spin parts of the Hamiltonian and the singlet- 
and triplet-superconductor responses are divergent. 

V. SUMMARY 

The one-dimensional 1/4-filled Hubbard model with alternating on-site interactions has 
been considered in the weak coupling approach using the rcnormalization group technique. 
The model is characterized by the decoupling of charge and spin degrees of freedom, but 
the dynamical nonequivalence of sites leads to the appearance of Umklapp processes in the 
system and to the dynamical generation of a gap in the charge excitation spectrum for 
Ua 7^ Ub, Ua > or Ub > 0. The ground-state phase diagram is obtained in the limit of 
second order renormalization. Depending on the sign and relative values of the bare coupling 
constants, there is a gap in the spin or charge excitation spectrum and the model system 
tends to superconducting or antiferromagnetic order at T = 0, with doubled period. It 
is shown by scaling to exactly soluble models that the terms corresponding to scattering 
processes, in which both incoming particles are from the same branch lead to the difference 
between the degrees of the divergency of CDW and SDW response functions. The role of 
interaction between particles on nearest and next-nearest neighbor sites is also considered. 
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TABLES 

TABLE I. The values of the fixed points for invariant couphng constants. The letters A...C 



describe the same regions as in Fig. 2 and 7 = (c + -^c^ + 8c) /2 



A B C D E 

gs{0) -2 -2 -2 

52(0) -1-7/2 1 1 

53(0) 2 2 -2 -2 

9p{0) 7 -2 -2 -2 -2 



TABLE II. The critical exponents ttj 

ABC D E 

Naib) -3/2+7/2 l(-3) 5 (-3/2) -3/2(5/2) -3(1) 

Xa(b) 5/2+7/2 1(5) -3/2(5/2) 5/2 (-3/2) 5(1) 

-3/2-7/2 1 5/2 5/2 1 

At 5/2-7/2 1 5/2 5/2 1 
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FIGURES 

FIG. 1. Classification of the scattering processes when interaction terms are expressed in a 

basis of right (sohd Hncs) and left (dashed hnes). a) Forward scattering in the vicinity of different 
Fermi points, b) Backward scattering, c) Umklapp processes. 

FIG. 2. The ground-state phase diagram of the model (1) {Vi = V2 = 0). The response 
functions corresponding to the phases shown in the parentheses have a lower degree of divergence 
than the others. 
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